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摘要 ”作为 聚变 反应 堆 候 选材 料 ， 低 活化 铁 素 体 钢 〈 马 氏 体 钢 ) 的 主要 成 分 是 Fe 和 Cr。 在 实际 工 况 下 ， 


聚变 反应 产生 的 氨 会 进入 材料 中 ,产生 微观 缺陷 ， 从 而 引起 辐 照 损伤 ， 影 响 材 料 的 抗 辐 照 能 力 ( 氨 脆 现 象 ) 


n 


及 其 它 物理 性 能 。 本 文通 过 第 一 性 原理 计算 ， 模 拟 得 出 反应 堆 主 体 材 料 Fe-Cr 合金 的 界面 ， 在 不 同 的 蔡 代 


位 和 间隙 位 迭 杂 He 原子 ， 对 其 结构 优化 后 计算 每 种 结构 的 形成 能 ， 得 出 He 原子 在 Fe-Cr 界面 不 同位 置 处 


NS 


的 稳定 性 。 与 此 同时 , 通过 对 体积 的 计算 得 出 He 在 靠近 界面 的 Fe 八 面体 间隙 位 时 会 引起 较 大 的 晶 格 畸变 ， 


根据 电子 态 密度 分 析 ， 推 断 是 He. Fe. Cr 三 种 原子 在 特定 位 置 相互 杂 化 影响 的 结果 ， 为 下 一 步 研究 核 级 


不 锈 钢 微观 缺陷 及 其 物理 化 学 性 质 提供 理论 依据 。 


关键 词 Fe-Cr 合金 ”第 一 性 原理 A 界面 


First principles study on helium doping at the interface of Fe Cr alloy 
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Abstract: As a candidate material for fusion reactors, the main components of low activation ferrite steel 
(martensitic steel) are Fe and Cr. In actual working conditions, helium produced by fusion reactions will enter the 
material, causing microscopic defects and radiation damage, affecting the material's resistance to radiation (helium 
embrittlement phenomenon) and other physical properties. This article simulates the interface of Fe Cr alloy, the 
main material of the reactor, through first principles calculations. He atoms are doped at different substitution and 
interstitial sites, and the formation energy of each structure is calculated after optimizing its structure. The stability 
of He atoms at different positions on the Fe Cr interface is obtained. At the same time, by calculating the volume, 
it was found that He would cause significant lattice distortion at the Fe octahedral gap near the interface. Based on 
the analysis of electronic density of states, it is inferred that this is the result of the hybridization of He, Fe, and Cr 
atoms at specific positions, provide a theoretical basis for the next research on the microscopic defects and 


physical and chemical properties of nuclear grade stainless steel. 
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根据 我 国 《 中 国 绿色 能 源 发 展 纲 要 〈2023-2033)》 和 “《“ 十 四 五 ”能 源 规划 》， 核 能 正 逐 


浙 成 为 我 国 重要 的 绿色 能 源 , 要 积极 有 序 发 展 核能 , 而 核反应 堆 材 料 对 于 氧 氮 
| 约 和 限 值 核能 发 展 的 因素 之 一 。 作 为 未 来 核 聚 变 反 应 堆 候选 材料 
一 低 活 化 铁 素 体 / 马 氏 体 钢 的 主要 成 分 就 是 Fe-9Cr 合金 "， 其 国内 外 主要 成 分 如 表 1 所 示 。 


其 一 直 是 | 
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的 抗 辐 照 能 力 
的 热门 材料 一 


而 在 工 况 下 ， 高 能 粒子 辐 照 会 让 材料 内 部 产生 缺陷 , 包括 空位 及 相关 原子 ， 


其 中 主要 包括 空 


位 ”、 间 隙 、 位 错 及 位 错 环 等 复杂 结构 ”"， 这 些 复 杂 结 构 会 破坏 材料 的 内 部 构造 ， 从 而 导致 


影响 材料 的 抗 辐 照 、 抗 蠕 变 等 能 力 。 因 此 , 本 文选 择 使 用 核反应 堆 中 最 常见 的 He 原子 作为 挫 


杂 原 子 ， 模 拟 其 在 Fe-Cr 合金 界面 时 的 情况 ， 从 而 进行 计算 分 析 。 
R1 国内 外 低 活 化 铁 素 体 / 马 氏 体 钢 的 主要 成 分 [4-9] 


Tab. 1 Main components of low activated ferrite / martensitic steel all around the world?! 


国家 型 号 主要 合金 元 素 微量 合金 元 素 
欧洲 EUROFER97 Fe (1%) ,9Cr W, V, Ta 
RS 9Cr2WVTa Fe 〈 基 体 ) ,9Cr 2W, V, Ta 
日 本 F82H Fe 〈 基 体 ) ,7.5Cr 2W，V，Ta 
日 美 联 合 JLF-1 Fe (GER) ,9Cr 2W, V, Ta 
日 美 联 合 JLF-2 Fe (GER) ,9Cr 2W, V, Ta 
日 美 联 合 JLF-3 Fe (4%) ,7Cr 2W, V, Ta 
中 国 CLAM Fe (基体) ,9Cr 1.5W, V, Ta 
中 国 CLF-1 Fe 〈 基 体 ) ,8Cr 1.5W, M, V, Ta 
中 国 SCRAM Fe (EK) ,9Cr 2W, V, Ni, Ti 
核 材 料 领 域 中 ， 研 究 人 员 针 对 于 Fe-Cr 合金 ， 主 要 通过 扫描 电镜 、 透 射电 镜 、 正 电子 淹 
没 等 方法 展开 了 微观 缺陷 研究 。 有 研究 表明 ，Fe-Cr 合金 中 Cr 原子 会 聚集 成 Cr 析出 相 "。 


Cr 原子 的 适量 添加 会 减弱 Fe-Cr 合金 中 位 错 及 空洞 的 形成 。 对 Fe-9%Cr 合金 进行 退火 实 


验 ， 当 温度 上 升 至 500C 时 ， 较 大 的 微观 缺陷 已 经 逐渐 消失 ， 当 温度 继续 升 


高 直至 800C 后 ， 


= 第 一 壁 


Cr 原子 开始 聚集 并 形成 Cr 析出 相 ， 最 终 溶解 。Fe-Cr 合金 作为 核反应 


P E 


的 基体 材 


料 , 研究 人 员 通 过 离子 注入 机 ， 电 化 学 等 多 种 方法 ,针对 其 微观 缺陷 内 部 的 氮 迁 移行 为 开展 
了 大 量 的 实验 研究 "9 ， 对 于 Fe-Cr 合金 中 的 所 行为 推断 入 总 结 如 下 ; 单个 或 多 个 氨 原 子 进 


入 材料 后 容易 与 空位 和 其 他 氨 原 子 相 结 合 ， 从 而 聚集 形成 氨 泡 ,这 一 现象 在 实验 中 
体现 出 来 ,同时 在 辐 照 的 作用 下 会 伴随 着 Frenkel 缺陷 对 产生 。 当 越 来 越 多 的 氨 原 子 进入 后 


能 很 好 的 


会 引起 品格 畸变 ,最 终 与 空位 不 断 吸引 形成 团 禾 ， 最 终 产 生肉 眼 可 见 的 氨 泡 ， 影响 材料 的 物 


理化 学 性 质 。 


本 次 研究 通过 采 月 


R 
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密度 泛 函 理论 ， 通 过 第 一 性 原理 计算 ， 可 以 高 精度 的 计算 得 出 He Ji 


子 在 Fe-Cr 界面 不 同位 置 的 稳定 性 ， 从 理论 上 得 出 He 原子 在 不 同位 置 时 对 Fe-Cr 界面 品格 
参数 的 影响 ， 为 核 级 不 锈 钢 界面 处 氨 迁 移行 为 的 研究 提供 理论 基础 。 


1. 计算 方法 


第 一 性 原理 又 名 从 头 计算 ,计算 特点 是 只 考虑 粒子 之 间 的 作用 原理 ， 通 过 量子 力学 的 方 
法 去 求解 薛 定 齐 方程 ,从 而 计算 得 到 材料 物理 化 学 性 质 的 方法 ,计算 过 程 中 不 含 任何 经 验 参 
数 。 本 次 使 用 第 一 性 原理 计算 软件 是 维也纳 从 头 计算 模拟 包 VASP) ”。 


NR 


2. 结构 优化 


本 次 优化 的 结构 模型 设置 为 Fe-Cr 界 面 ， 属 于 典型 的 体 心 立方 结构 。 对 于 和 氨 原 子 返 杂 的 
位 置 设 置 为 蔡 代 位 和 间隙 位 ,其 中 间隙 位 又 分 为 八 面体 间隙 和 四 面体 间隙 。 本 论文 使 用 VASP 
软件 对 Fe-Cr 界 面 进行 结构 优化 ， 优 化 方法 采用 6GA《〈 交 换 关 联 泛 函 ) 和 PAW (级 加 平面 波 ) 
下 。 另 外 设置 ISIF=3，ENCUT=400 eV， 直 到 优化 至 每 个 原子 受 力 二 0. 001 eV。 优 化 结束 后 


继续 进行 静态 计算 ，IBRION=-1， 得 出 界面 处 原子 分 布 模型 (图 1)。 


1 Fe-Cr 合金 界面 的 原子 模型 (1,2,3,4 分 别 代表 平 行 于 xoy 平面 的 切面 ，Fe 原子 用 黄色 表示 ，Cr 原子 用 


蓝 色 表示 ) 
Fig. 1 Atomic model of Fe-Cr alloy interface (1, 2, 3, 4 represent the sections parallel to the xoy plane, respectively, 


yellow represents the Fe atom and blue represents the Cr atom.) 


S 


中 ，Cr 原子 用 蓝 色 表示 ，Fe 原子 用 黄色 表示 ， 切 线 1、2、3、4 分 别 表 示 He 原子 


所 在 的 平面 。 另 外 ，T (Tetrahedron) 表示 He 原子 位 于 四 面体 间隙 位 ，0 COctahedron) 


表示 He 原子 位 于 八 面体 间隙 位 、S CSupersede) 表示 He 原子 位 于 替代 位 ， 为 了 进行 对 


比 ， 另 加 入 空位 : V (Vacancy) 表示 空位 在 蔡 代 位 (图 中 某 个 Fe 或 Cr 原子 消失 ， 其 原 


所 在 的 位 置 没 有 任何 原子 )。 比 如 : T3 就 表示 He 原子 位 于 3 号 切面 上 ， 并 处 于 由 一 个 Cr 
原子 和 三 个 Fe 原子 所 构成 的 四 面体 的 中 心 。 在 结构 优化 的 基础 上 ， 分 别 把 He 原子 添加 
到 以 上 位 置 处 ， 再 对 其 进行 结构 优化 ， 优 化 参数 与 之 前 相同 ， 优 化 后 再 进行 一 次 静态 计 
算 ， 并 利用 以 下 公式 计算 其 形成 能 ”: 


Erormation 一 Ere. Cr. He 一 Eg. Cr 一 EHe (1) 


E Formation = Ere. Cr.He 一 Ep. Cr 一 EHe + Esingle Fe/Cr (2) 


EB. Er o n RESP ZRUH. He 原子 的 Fe-Cr 界面 体系 的 总 能 量 ，E;.. 6 表示 Fe-Cr 界 


E. 
里 
面体 系 的 总 能 量 ，B 表示 单个 He 原子 在 一 个 相对 无 穷 大 空间 内 的 总 能 量 ”。Eswsuwe 表 


示 被 蔡 代 的 单个 Fe 或 Cr 原子 的 静止 能 量 。 在 经 过 相关 计算 过 后 ， 得 
界面 各 个 位 置 的 形成 能 如 图 2 所 示 。 


来 的 He 在 Fe-Cr 


Er 


形成 能 (eV) 


四 面体 间隙 
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FeCrV1 V2 V3 V4 S1 S2 S3 S4 O1 O2 O3 O4 T1 T2 T3 T4 - 
He 和 空位 所 在 的 不 同位 置 


图 2 He 和 空位 在 Fe-Cr 合金 界面 不 同位 置 处 的 形成 能 


Fig. 2 Formation Energy of He and Vacancies at Diffe rent Positions of Fe-Cr Alloy Interface 


3. 形成 能 


从 图 中 y1，V2，V3 及 V4 的 形成 能 可 以 得 出 ， 当 Fe-Cr 界面 中 少 了 一 个 原子 ， 导 致 
单 空位 产生 时 ， 其 结构 是 不 稳定 的 。 相 比 于 Cr 原子 处 出 现 空位 ， 空 位 更 倾向 于 在 远离 界 
HR Fe 原子 位 置 形成 V4 位 置 表示 远离 界面 的 Fe 原子 位 置 ， 其 形成 能 相对 最 小 ， 稳 定 
性 也 相对 最 高 )。 换 句 话 来 讲 ， 当 Fe-Cr 界面 中 出 现 单个 空位 时 ， 空 位 更 倾向 于 在 远离 界 
看 的 Fe 原子 位 置 处 形成 ， 即 四 周 都 是 Fe 原子 的 位 置 。 这 为 接 下 来 研究 Fe 基 合金 核反应 
堆 材料 中 空位 的 形成 位 置 ， 提 供 了 理论 依据 。 

根据 图 中 He 原子 在 不 同位 置 ($S1，S2，S3，S4，01，02，03，04，T1，T2，T3，T4) 
的 形成 能 可 以 得 出 ， 氨 原子 在 四 面体 间隙 位 时 相对 最 稳定 ， 其 次 是 在 八 面体 间隙 位 ， 在 
蔡 代位 时 相对 最 不 稳定 。 在 这 之 中 ， 氨 原子 在 T1 位 置 处 的 形成 能 最 低 ， 在 S1 位 置 处 的 
形成 能 最 高 。 由 此 可 以 得 出 ， 当 单一 He 原子 注入 Fe-Cr 界面 后 ,会 首先 倾向 位 于 远离 界 
看 的 Cr 四 面体 间隙 位 中 ， 如 有 多 个 He 原子 进入 ， 会 进入 其 他 四 面体 间隙 位 ， 直 至 靠近 
界面 的 Fe 四 面体 间隙 位 都 填 满 后 ， 才 会 进入 八 面体 间隙 位 中 ， 一 直到 八 面体 间隙 位 也 填 
满 ， 过 多 的 He 原子 就 会 开始 替代 原 有 的 Cr 原子 与 Fe 原子 ， 取 代 他 们 的 位 置 ， 形 成 氨 泡 
或 氮 空 位 复合 体 ， 从 而 产生 微观 缺陷 。 对 于 被 蔡 代 的 Fe 原子 和 He 原子 ， 由 于 脱离 了 原 
有 的 唱 格 位 置 ， 会 因为 较 大 的 原子 半径 ， 对 原 有 晶 格 参数 产生 较 大 影响 ， 在 材料 中 形成 
微观 缺陷 ， 比 如 空位 、 空 位 团 、 团 艇 、 缺 陷 对 、 位 错 、 位 错 环 等 ， 进 一 步 影 响 材 料 的 物 
理化 学 性 质 。 


4. 体积 变化 


为 了 更 直观 的 得 出 He 在 Fe-Cr 界面 不 同位 置 的 影响 ， 第 一 性 原理 计算 得 出 He 和 空 
位 在 Fe-Cr 界面 不 同位 置 处 的 体积 如 图 3 所 示 。 
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图 3 He 和 空位 在 Fe-Cr 合金 界面 不 同位 置 处 的 体积 


Fig. 3 Volume of He and Vacancies at Different Positions of Fe-Cr Alloy Interface 


显而易见 ， 相 比 于 退火 后 没有 微观 缺陷 的 Fe-Cr 合金 界面 , 无 论 是 He 在 替代 位 还 是 
空位 在 蔡 代 位 ， 都 会 使 得 Fe-Cr 界面 的 体积 相对 减 小 ， 且 单 空位 的 出 现 ， 更 会 明显 降低 
Fe-Cr 界面 的 体积 。 当 He 原子 出 现在 间隙 位 时 ， 会 明显 的 增加 Fe-Cr 界面 的 体积 。 值 得 
注意 的 是 ， 当 He 原子 出 现在 Fe-Cr 界面 的 不 同位 置 时 (01，02,，03，04,， T1, T2, T3, 
T4), T3 位 置 的 体积 比 其 他 位 置 的 体积 都 大 ， 且 体积 的 变化 近乎 是 别 的 位 置 的 两 倍 。 这 

计算 结果 表明 ， 当 He 原子 位 于 靠近 界面 的 Fe 八 面 体 间 隐 位 时 ， 会 引起 整个 界面 出 现 
较 大 的 唱 格 畸变 ， 这 种 畸变 会 一 定 程度 上 改变 材料 的 形状 ， 从 而 影响 材料 的 物理 化 学 性 
质 。 


ed 


针对 于 这 一 位 置 的 特殊 现象 , 我 们 继续 对 He 在 此 位 置 的 结构 模型 进行 了 多 次 精度 更 
高 的 优化 ， 其 优化 结果 都 在 358. 94 —359. 2A ZR] (1A=10 "m) 。 与 此 同时 ， 使 用 p4vasp 
软件 得 出 了 它 的 态 密度 图 以 及 包含 不 同 电子 轨道 的 分 波 态 密度 图 ， 并 与 其 他 八 面体 间 辽 
位 《01, 02, 04). 的 电子 态 密度 图 进行 了 比较 《图 4)。 


Density (states/eV) FeCrHe 


-10.0 -8.0 -6.0 -40 -20 0.0 2.0 
Energy (eV) 


Density (states/eV) FeCrHe 


60.0 
^| He 在 02 位 置 
50.0 e 
45.0 
40.0 
35.0 
30.0 
25.0 
20.0 
15.0 
10.0 
EA 
0.0 
-20.0 -18.0 -16.0 -14.0 -12.0 -10.0 -8.0 -6.0 -4.0 -2.0 0.0 2.0 
Energy (eV) 


Electronic Density of States 
Density (states/eV) FeCrHe 


—————hhNJhJOhAhÓPRACO C9 CO 


ONPOoOONPROOONPROOONR 
OOOOoOoOoooooooooooo 


Density (states/eV) FeCrHe 


He 在 04 位 置 


-20.0 -18.0 -16.0 -14.0 -12.0 -10.0 -8.0 -6.0 -4.0 -2.0 0.0 2.0 
Energy (eV) 


xm], 2À 


WR] 
o 
> 
- 
En 
3 
E 
an 
[NY 
时 


图 4 He 在 Fe-Cr 合金 界面 不 同 八 面体 间隙 位 01，02，03，04 的 电子 态 密度 


坐标 表示 电子 态 密度 


(03 中 红 绿 蓝 三 条 曲线 分 别 表 示 Fe、Cr、He 原子 的 分 波 电子 态 密 度 ) 


Fig. 4 Electronic Density of State Charts of He at Different Octahedral Gap Points O1, O2, O3, O4 at Fe-Cr Alloy 
Interface (The three curves of red, green and blue in O3 represent the fractional density of states of Fe, Cr and He 


atoms, respectively ) 


电子 态 密 度 能 反应 得 到 电子 在 特定 能 量 区 域内 的 分 布 情况 ， 其 态 密 度 图 能 反应 电子 
在 不 同 轨 道 的 分 布 ， 从 而 反映 出 原子 之 前 的 相互 作用 。 从 图 4 中 可 以 看 出 ， 其 他 八 面体 
间隙 位 〈01，02，04) 在 能 量 -20Ev~-18eV 处 有 一 个 波峰 ， 而 He 在 靠近 界面 的 Fe 八 面 


体 间隙 位 C03) 时 的 波峰 出 现在 能 量 -18eV~-16eV 之 间 。 波 峰 的 偏 移 表 示 He 原子 在 此 位 
置 会 与 周围 的 Cr 原子 和 Fe 原子 有 特殊 的 相互 作用 , 该 峰 位 不 只 是 一 个 原子 作用 的 结果 ， 
而 是 Fe、Cr、He 三 种 原子 所 代表 的 红 绿 蓝 三 条 线 都 有 波峰 产生 ， 说 明 He 原子 的 电子 轨 
道 与 Fe 原子 和 Cr 原子 的 电子 轨道 出 现 了 杂 化 现象 。 该 峰 是 由 Fe, Cr, He 所 共同 组 成 的 
杂 化 峰 ， 是 Fe. Cr. He 三 原子 之 间 一 种 特殊 的 相互 作用 。 但 是 相 比 于 其 他 八 面体 间隙 位 
来 说 不 稳定 ，He 原子 很 容易 迁移 到 别 的 八 面体 间隙 位 以 及 四 面体 间隙 位 。 


第 一 性 原理 计算 得 出 ， 当 He 原子 进入 Fe-Cr 合金 界面 时 ， 首 先 倾向 于 向 远离 界面 位 
置 的 Cr 四 面体 间隙 位 置 处 移动 , 当 这 一 位 置 被 完全 占据 , 之 后 依次 是 其 他 四 面体 间隙 位 ， 


八 面体 间 隐 位 ， 蔡 代位 。 这 一 结论 有 助 于 进 
而 研究 材料 的 抗 辐 照 能 力 。 

当 He 原子 进入 靠近 Fe-Cr 合金 界面 的 Fe 八 面体 间隙 位 时 , 会 引起 较 大 的 唱 格 畸变 ， 
这 是 He、Fe、Cr 三 种 原子 在 特定 位 置 相互 杂 化 影响 的 结果 。 这 使 得 界面 的 原子 排列 发 4 
较 大 改变 , 从 而 影响 材料 的 微观 缺陷 结构 , 该 发 现 为 下 一 步 研 究 Fe-Cr 界面 内 He 在 Fe 八 
用 体 间 隙 位 时 对 合金 造成 结构 影响 提供 理论 依据 和 研究 思路 。 


SE 


研究 核反应 堆 材料 中 的 氨 迁 移行 为 ， 从 
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